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Far-infrared bolometric detectors are used extensively in ground-based and 
space-borne astronomy, and thus it is important to understand their optical 
behaviour precisely. We have studied the intensity and polarisation response of 
free-space bolometers, and shown that when the size of the absorber is reduced 
below a wavelength, the response changes from being that of a classical optical 
detector to that of a few-mode antenna. We have calculated the modal content 
of the reception patterns, and found that for any volumetric detector having a 
side length of less than a wavelength, three magnetic and three electric dipoles 
characterise the behaviour. The size of the absorber merely determines the 
relative strengths of the contributions. The same formalism can be applied to 
thin-film absorbers, where the induced current is forced to flow in a plane. In 
this case, one magnetic and two electric dipoles characterise the behaviour. 

The ability to model easily the intensity, polarisation, and straylight charac- 
teristics of electrically-small detectors will be of great value when designing 
high-performance polarimetric imaging arrays. 
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1. Introduction 

Focal-plane arrays of far-infrared bolometers are used extensively in ground-based and space- 
borne astronomy. They can be divided into two categories depending on how power in the 
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incident radiation field is delivered to the sensing element. In horn or antenna-coupled bolo- 
meters, power is collected by a metallic antenna and then delivered to the sensing element 
through a transmission structure. This transmission structure could be either a single-mode 
transmission line terminating in a resistive load on the bolometer [1], or a waveguide in 
which the bolometer is mounted [2]. In free-space coupled bolometers, power is absorbed 
directly in a bulk material or resistive film [3-5] . The optical behaviour of single-mode horn 
and antenna-coupled bolometers is well understood. The antenna couples power from the 
multiple modes of the incoming partially coherent radiation field into the single mode as- 
sociated with the transmission structure that connects to the antenna. The fully coherent 
antenna pattern can be referenced through the optical system onto the sky, to give the field 
distribution to which the complete instrument is sensitive. The equivalent analysis for mul- 
timode horns (such as Winston cones) and free-space bolometers is more subtle because the 
absorber is receptive to power arriving in a number of optical modes simultaneously. In fact, 
the natural modes of the optical system, the natural modes of the field travelling through 
the optical system, and the natural modes of the detector are all distinct, and it is more 
complicated to calculate the precise form of the field distributions to which an instrument 
responds. 

In order to optimize the design of astronomical instruments it is vital to understand 
the behaviour of free-space bolometers. An optical theory of bolometric detectors has been 
proposed by Withington [6] and Saklatvala [7]. According to this theory, bolometers are 
characterised with respect to the field on a reference plane, which is deemed to be the 
collecting surface of the device. In the case of antenna-coupled bolometers, the collecting 
surface may be the aperture of a horn, whereas in the case of free-space bolometers, it may 
be the surface of a resistive film. It is then possible to expand the incident field over this 
surface in terms of a set of vector basis functions. In the case of partially coherent fields, the 
expansion coefficients are random variables. For a particular source distribution, however, it 
is possible to find a set where there is no correlation between the expansion coefficients. We 
shall refer to the members of this set as the natural modes of the incident field. Similarly, it 
is possible to characterise the response of a bolometer in terms of a set of natural modes, 
meaning that each mode contributes independently to the recorded output. The total power 
measured depends on how the natural modes of the field scatter into the natural modes of 
the detector. 

A critical question is how do we determine the optical modes to which a system of absorbing 
structures is sensitive? The approach often taken at optical wavelengths is to treat a detector 
as being sensitive to the intensity of the field integrated over the area of the absorber. 
This approach corresponds to assuming that the response is fully-spatially incoherent and 
unpolarized over the physical extent of the absorbing surface. To determine the reception 
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pattern on the sky, a geometrical optics approach is then taken, meaning that the system is 
assumed to be unpolarized and sensitive to intensity only over the re-imaged outline of the 
pixel. In this paper we show that this approximation is only valid when the wavelength A 
is much smaller then the physical size p of the absorber. This condition is usually satisfied 
by optical instruments, because the typical size of a CCD pixel is 5 x 10 _6 m compared 
with 5 x 10 _7 m for the wavelength of visible light: p/X > 10. At far-infrared and sub-mm 
wavelengths, however, devices are being constructed with absorbers that are smaller than 
the wavelength of operation: p/X < 1. In order to accommodate this more general case, a full 
theoretical model of power absorption has been developed [8] , which takes into account the 
state of coherence of the incoming field, and the state of coherence of the physical processes 
responsible for absorption. The scheme is related to the Method of Moments [9]. Although 
the approach is rigorous, it is too cumbersome for certain design problems, particularly when 
one wants to use physical optics to model the behaviour of large-format polarimetric arrays. 

In this paper, we adopt an intermediate position, and describe the optical behaviour 
of free-space bolometers in terms of modes. The underlying principles have already been 
presented [10, 11], but here we extend the analysis considerably. The basic idea is to place 
a bolometer at the back of an ideal imaging system, which characterises the fact that an 
incoming field has to travel through free space before arriving at the detector. The degrees 
of freedom associated with a propagating field are limited, and these restrictions introduce 
spatial correlations and cross polarisation into the reception pattern. Indeed, as the pixel- 
size is reduced below a wavelength, the reception pattern becomes increasingly polarized and 
divergent from its re-imaged geometrical shape. A key part of the work described here is that 
we explore the state of coherence of the angular reception pattern of a naked pixel. In this 
way, we show that the behaviour of small volumetric detectors, p/X < 1, can be characterised 
by three magnetic and three electric dipoles whose relative strengths are determined by the 
size of the absorber. In the case of thin-film absorbers, where the induced currents are forced 
to flow in a plane, only one magnetic and two electric dipoles are present. 

2. Bolometric detection 

The problem of calculating the power absorbed by a free-space bolometer from a field in any 
state of spatial coherence has been addressed by Withington [6] and Saklatvala [7] . Assume 
that the incident field is temporally stationary, so that different frequency components, u ;, 
can be treated independently. Let T> be the input reference plane of the detector, and let 
r t be a position vector on this plane. Let Ep(r t ,w) be the analytic-signal representation of 
the incident electric field over V. The first-order spatial correlations can be described by the 
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space-domain correlation dyadic 


Sv(r n ,r n ,u)= (Ei)(r tl ,w) ~E* v (r t2 ,u)), (1) 

where { ) denotes the ensemble average. £(r tl ,r tl ,u) has physical significance beyond the 
coordinate system used, but a matrix representation can be found by projecting £(r tl , r^, u) 
onto a specific basis: £ij(r n , r n , u) - x z (r tl ) ■ £ v (r tu r t2 ,u) ■ itj(r t2 ). 

Withington and Saklatvala [6, 7] have shown that the power absorbed by any bolometer 
can be written in the form 

(p(w)) = / d 2 r tl / d 2 r t2 Dl(r n ,r t2 ,u) ■ ■ £ v (r n ,r t2 ,u), (2) 

Jv Jv 

where the double-dot denotes tensor contraction to a scalar, and f denotes the adjoint. In (2) 
a new dyadic field, Dj)(r a ,r t2 ,u), has been introduced, which characterises the relationship 
between the state of coherence of the incident field and the output of the detector. We shall 
refer to it as the detector response function. 

Both £v( r tu r t2, and T> D (r t i,r t2 ,uj) are Hermitian in all variables, and admit diagonal 
decompositions: 


&(r,i,r H .w)= U‘“»(r K ,w). (3) 

n 


{A^} is a set of scalars, and {U^ n ^(r t , a;)} a set of vector fields. Let {A^} and {UH r t> w )} 
be associated with the decomposition of £ D (r t \,r t2 ,u;), and H" 1 )} and {R( m )( rij a;)} with 
the decomposition of V D (r a , r t2 , u>). Substituting these representations into (2) gives 


(pH) 


J> (m) A (n) I I d 2 r t R/ m >*(r t ,u) ■ U^(r t ,u) 

m,n ^ 


( 4 ) 


(3) corresponds to finding the coherent modes of the field [12], with the A^ corresponding 
to the fraction of the total power carried by each mode. Each Rd m )(rj , u ) is a fully coherent 
field to which the detector is sensitive, and each 7^ is the effective responsivity of the mode. 
In the case of antenna-coupled bolometers, only one 7^ is non-zero, and (4) takes on the 
well-known form for the power coupled between a coherent field and a radio antenna. 


3. Bolometric imaging 

Consider an optical system comprising two ideal thin lenses separated by the sum of their 
focal lengths: Figure 1. A planar bolometer is placed at the focal- plane, T , of one of the 
lenses. The system could represent an ideal astronomical telescope, with suitable angular 
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scaling at the source plane, S, or a microscope of some kind. The purpose of the model is 
to bring out fundamental restrictions on the imaging properties of bolometers, aside from 
constraints due to the optical system. We therefore assume that the lenses are large compared 
with the wavelength, and that they do not introduce cross-polar scattering. Obviously, the 
system cannot image the whole field at S, only the component that is able to propagate to 
T. 

Assume that the response function of the bolometer at T is known, Djr(r tl: r t2 ,uj), but 
we would like to characterize the instrument’s imaging capabilities at <S. Therefore, it is 
necessary to relate the response functions at T and S. Prom (2), the power absorbed is 

<p(w)> = J <?r n J d?r t2 D^{r n ,Y t2 ,u) ■ ■ S v (r tl ,r t2 ,uj), (5) 

where £p(r tl , r t2 , u>) is the coherence tensor of the propagating part, E P (r t , u), of the source- 
plane field, E,$(rt,a;). The propagating part can be found by using the angular spectrum of 
wavefields [13]. If the Fourier transform of Es(r t ,u>) is 

A(k t , u) = ~ J d 2 r t E 5 (r*, u) exp(— • r t ), (6) 

then each component is associated with an expansion function of the form A(k t , u>) exp(zk t • 
r t ) exp(ik z (li t )z), where k = k t + k z z, and |k t | 2 + (k z (k t )) 2 — ( u/c ) 2 ; it is assumed the 
source plane is at z = 0. Because we are only interested in the part of the field that is 
able to propagate from S to J 7 , we must consider only those Fourier components for which 
|k t | < (u/c), and those having vector components transverse to the direction of propagation. 
We therefore simply discard those Fourier components not associated with plane waves, and 
invert the transform. Strictly speaking, this approach ignores certain boundary conditions 
on the tangential component of the A-field at S, but it is appropriate for our purposes. 
Applying these constraints gives 

E P (r t , u) = J d 2 \t t (l - kkj • A(k t , u) exp(ik t • r t ), (7) 

\*A<t 

where T is the unit dyadic in three-dimensions, and X — kk strips out any longitudinally 
polarized component. Using (7), the correlation dyadic of the propagating field becomes 

! p (: rti,ra,w) = J J d ^ t2 ex P(* k « • r n) exp(-ik t2 ■ r t2 ) 

|k*l<f l k *l<“ (8) 

x (l — kiki^ • A(k tl ,k t 2,uj) ■ (x — k 2 k 2 ^ , 
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where *4(k tl , k t2 , to) — { A(k tl , uj)A*(k t 2 , u>) ) is the ^-domain correlation dyadic. 

It is now possible to relate the detector response functions at the focal and source planes. 
Substituting (8) in (5), and changing the order of integration, gives 

(p(w)) = Jd 2 k t i Jd 2 k t2 (t — kiki j • ^(kti,k t 2 ,u;) • (x - k 2 k 2 )] f • • „I(k u ,k t 2 ,u;). 




( 9 ) 

^(k t i,k t2 ,a;) is the /c-domain representation of the detector response function on the focal 
plane: 

^(kti,ka,w) = — ^ f Sr n f d 2 r t2 /^(r^r^w) exp(-ik tl -r tl ) exp(+ik t2 -r t2 ). (10) 

Finally, (9) can be cast back into the space domain giving 

(p(w)) = / / d 2 r t2 Dl(r n ,r t2 ,uj) ■ ■ £s(r tl , r t2 , w), (11) 

./S 


with 


£>s(rti,r t 2 ,u;) = J d 2 k n J d 2 k t2 (l - kiki) • £?(k n ,k t 2 ,a;) • (x-k 2 k 2 ) 


|kt|<“ |k»l<* 

x exp(ik« • r tt ) exp(-ik t2 • r t2 ). 


( 12 ) 


(10) and (12) together provide a relationship between the response functions at T and 
S. Notice that the state of polarization and spatial-coherence of the detector’s response 
are altered when referenced to S. The limits on the integrals filter the spatial-frequency 
components, and introduce coherence into the reception pattern, which affects its modal 
structure. Although an optical system will, in reality, introduce additional constraints, (12) 
represents the minimum filtering that a re-referenced response function must incur to give a 
reception pattern. 

In astronomy, sources are incoherent, and the equivalent source-plane field can be written 

2 

£s(r n ,r t2 ,Lj) = ^ E ij ( r *i da^) x*Xj 5 (r« - r t2 ) , (13) 

i,j= 1 

where {xj} is a set of orthonormal basis vectors lying in the plane. Substituting (13) into 
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(11) and letting — 5q • D s ■ x^-, we obtain 

2 r 

(p(w)) = / d 2 r t D^rt.rt.u) E i:j (r t ,r t ,u). (14) 

i,3 = 1 

It follows from (1) that the Eij(r t , r t , u) are elements of a 2 x 2 Hermitian matrix. The matrix 
can therefore be expanded in Pauli-spin matrices, 


{ 


< 7 ^ 





or 

E i:j ( r t ,r t ,u) 

where the coefficients are given by 


\ Ex“W)4’. 


X {k) (r*, u}) = ^2 Eij{ r t , r t , u) a$ 


(*) 


*j'=i 


(15) 


(16) 


(17) 


If we let {xi, x 2 } = {x, y}, then the correspond to the Stokes parameters of the field on 
the source plane. Substituting (16) into (14), we obtain 

(p(w)) - J d 2 rt 1V(r t , u) • xfa, u), (18) 

where x is the vector of Stokes parameters. The components of r t , uj) are given by 

Ti {k) (r u ijj) = \Y1 r t, UJ ) <T ji\ (19) 

*J=1 

and they are the position-dependent Stoke’s parameters of the partially coherent reception 
pattern of the bolometer. Thus, in the case of spatially incoherent fields, the full response dy- 
adic can be replaced with a position-dependent vector, x( r t, w), which describes the response 
of the detector to different states of polarisation. 


4. Square bolometers 

Consider a square bolometer having side length p, and assume that the physical process 
responsible for energy absorption is uniform, unpolarized and local. The focal-plane response 
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function is then 


XV( *,.*,«) = ( at5 ^ ~^ W- W- W-W * 1 . (20) 

I 0 elsewhere 

r t = (x, y ), and we have assumed that the origin lies at the centre of the absorber. X and 0 
are the 3x3 identity and zero dyadics. a is a model parameter representing the absorptivity 
of the bolometer and is assumed to be a constant. At T, the natural modes of the bolometer, 
ignoring edge effects, form a degenerate continuum spanning all possible field states over the 
region defined by the absorber. An important question is whether this model can represent 
any physically realizable device. The model was formulated originally to describe black-body 
absorbers, and so can best be thought of as representing the response of a square aperture in 
the wall of an integrating cavity. Alternatively, it can be thought of as representing a thick, 
volumetric absorber, where scattering has, strictly speaking, been ignored. The absorber is 
volumetric in the sense that it responds to the ^-directed component of the incident electric 
field, or equivalently to current flow in the z direction, as well as in the plane. By considering 
ohmic power dissipation in the absorber, it is very easy to arrive at (20). This model contrasts 
with that of a thin film absorber, which will be considered later. 

Combining (20), (10) and (12), gives the response function on the sky, 

t>s(r a ,r a ,ui) = J 'I'ku J d 2 k, 2 M (kn.ka.w) exp (-Sk tl ■ r a ) exp (+ik 12 ■ r t2 ) , 

( 21 ) 

where 

^(k«,k t2 ,u;)= |k tl |) H(“- |kfl|) 

xc^J-kiki) • (x - k 2 k 2 ) (22) 

x S(k n ,k t2 ,u) , 

and H is the Heaviside step function. The two Heaviside terms represent the necessary 
restriction on |k t | for k z to be real and positive. The dyadic operators in (22) strip out 
any longitudinal field component, and introduce coupling between the polarisation states of 
the wave-fields. We shall refer to S (k t i,k t2 ,u;) as the shape function because it encodes all 
information about the shape of the pixel. It is given by 

S (k n , k*a, u ;) = — ^ f d 2 r t exp(i(k tl - k t2 ) • r t ) (23) 

(27T) Ja 

where A is the square region occupied by the absorber, |x|, \y\ < p/2. Tree-space absorbers 
having different shapes can be modelled simply by changing A. In the case where A is a 
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square 


S (k u ,k t 2 ,a;) = 


/ 


cZx exp (i(k n - k* 2 ) • x x) 


x 


(2tt)2 

' r\ 

/ ^ dy exp(*(k tl -k t2 )-yy) 


= ( l ?) 2 sinc [^ (k,2 " k,,) ' ii ] 

x sine (kt 2 - k ( i) ■ yj , 


( 24 ) 


where we have introduced the normalised sine function sinc(x) = sin(7rx)/(7rx). 

Although (24) has a simple form, there is no obvious analytic solution to (21), and so the 
transformation must be achieved numerically. There are two possible options: The first is to 
evaluate the integrals directly. In this case, the symmetries can be exploited, meaning that 
only the values of the coherence dyadic over an octant of the plane about the pixel centre 
need to be calculated explicitly. The second, which we will describe here, involves calculating 
the response function using a superposition of fundamental modes. It will be shown that 
this latter technique is both numerically advantageous and provides better physical insight, 
especially in the p/X limit in which the number of required modes is small. 

Prom (22) and (24) it can be seen that M. (k w , k t2 , u) is Hermitian in all variables, and 
therefore has a decomposition of the form 


•M(k tl ,k t2 ,w) =^7 (n) V (n >(k tl ,w) Vl*l*(k a ,u), (25) 

n 

where the 7 ^ and (r t ,u) are solutions of 

J d 2 k t2 M(k tl ,k t2 ,u) . y( n ) (k t 2 ,u>) — 7 (n) V<">(k n,u). (26) 

The integrals are evaluated over the entire k r plane. The (r ti , u) are the natural modes 
of the detector, as discussed in Section 2, but now in the fc-domain, and for fields on the sky. 
Substituting (25) into (21), gives 

T> s ( r t i,r t2 ,u) ^ 7 W U W (r t , W )UW*(r ilW ). (27) 

71 

where 

U (n) (r t , uj) = ^ J d 2 k t (k*, u) exp (ik t ■ r t ) . (28) 
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It follows from (27) that the 7^ and (r t , u) are solutions of 

J d?r t V s (r tl , r t2 ,u) • U (n) (r i2 , u) = 7 (n) U (n) (r n , u) , (29) 

which are the natural optical modes of the detector’s response function on the source plane. 
This analysis shows that the modes constitute a set of domain-independent fields, and there- 
fore have true physical significance. For our purposes, it provides an efficient alternative to 
(21) for transforming between M. (k tl , k t2 ,a;) and V s (r tl , r t2 , u). The procedure is as fol- 
lows: (i) Calculate a sampled version of M. (k tl , k t2 , a;) using (22); (ii) diagonalise according 
to (26), (iii) transform the modes into the space domain using (28); and (iv) reconstruct 
V s (r n ,r t2 ,u) using (27). 

This modal approach has two distinct advantages over the direct evaluation of the integral 
in (21). First, though the modal decomposition requires some numerical overhead, under- 
standing the natural modes of the system provides significant physical insight into its optical 
behaviour, especially in the limit where the pixel is smaller then the wavelength. Second, 
the modal approach is numerically more efficient. Consider transforming the response func- 
tion from a grid of points in the Fourier-domain to a grid of points in the space-domain. 
The number of two-dimensional Fourier transforms that must be performed to evaluate (21) 
scales with the square of the number of sample points. In the modal case it scales with the 
number of modes with nonzero eigenvalues. Even when pixel symmetries are exploited, in 
most cases the modal approach requires fewer transforms then the direct approach. 


5. Numerical implementation 

To implement the scheme it is necessary to sample Ai(k tl , k t2 , lj). Consider sampling 
vVl(k t i, k t2 , oj) over a set of points jk^ j forming a square-grid over part of the k* plane. 
Assume that the grid is centred on the origin, and consists of N points in total, spaced 
equally in the k x - and ^-directions at intervals of 27 t/L: Figure 2a. (26) then becomes a set 
of simultaneous linear equations: 



X, ■ M( kj r) ,kj' ) ,w) • x„ x»-V < ’ ,, (ki' , ,w) 



(30) 


for q — 1 to 3 and r = 1 to N. (2i t/L) 2 is the area of the plane associated with each sample 
point, and {x 9 } is any set of basis vectors spanning three-dimensional space. Let the matrix 
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representation of (30) be 


3 N 

y^f ) = 7 ( "V 1 ( ” , , (31) 

3 = 1 

then the problem of solving (26) is reduced to finding the eigenvalues, 7 ^ n \ and normalised 
eigenvectors, , of the matrix M. Computationally, we calculate the M ^ using (22), solve 
for the eigenvalues and eigenvectors using a library function, and then extract the modes 
sampled over the grid. Because of the choice of grid of points, the Fourier transform reduces 
to a Discrete Fourier Transform (DFT), for which optimised routines can be used. The space- 
domain response function, sampled over the grid of points shown in Figure 2b, can then be 
reconstructed from the transformed modes using (27). 

When finding the natural-modes numerically, we must sample the k 4 plane sufficiently well 
to account for all degrees of freedom. It is easy to appreciate that unless we increase the total 
number of points, N, as we increase the sampling rate in the fc-domain, the resolution in the 
space domain decreases. In fact, to maintain resolution as the ^-domain spacing is decreased, 
N must scale as L 2 . Increasing N, however, comes with the significant numerical overhead 
of finding the eigenvector decomposition of a large matrix. The problem can be mitigated 
to some extent by remembering that all the elements of M in rows or columns where either 
k a > w/c or k tt > u/c are zero. It is straightforward to show that these blocks of zeroes 
do not affect the eigenvector decomposition, and further that for all modes V^(k t ) — 0 
when k t > u/c. Consequently, it is possible when performing the eigenvector decomposition 
to consider the reduced matrix M s consisting only of pairs of points where k tl < u/c and 
k *2 < oj / c. The values outside the domain k t can then be filled with zeroes in a later step to 
increase the resolution in the spatial domain when the mode is Fourier transformed. 

A consequence of this result is that the maximum number of modes that can be recovered 
depends only on the row (or column) length of M s , which is equal to three times the number 
of points in the grid that fall inside |k t | < u/c, and is increased by increasing L. For the 
discretization to be a good approximation to (26) the number of degrees of freedom must, 
at the very least, exceed the actual number of natural modes of the system. In practice, this 
number is not known in advance. To ensure all natural modes of the system were accounted 
for, we iteratively ran the code with increasing values of L for each value of p/X. In each 
case, it was assumed L was sufficient when the non-zero eigenvalues stopped increasing in 
number and converged on final, individual, values. 

6. Simulations I 

Figure 3 shows the components of lZ(r t ) for p/X — 0.25, 0.5, 1, 2 and 4. The black square in 
the centre of each plot represents the geometrically re-imaged absorber. For each value of 
p/X , the plots have been normalised to the peak value of TT-i, and 7Z\ is plotted against a 
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single scale spanning 0 to 1 for all p/X. The plots for 7Z 2 to 7Z^ have been plotted using a 
common colour scale for any given p/X, but the scale changes for different values of p/X. 
In [11], Chuss et. al. presented the same quantities using a different, less efficient, numerical 
method, and it is pleasing to see that the results are in exact agreement. 

The IZi characterise the response of the bolometer to spatially-incoherent partially polar- 
ized fields. IZi represents the responsivity to the Stokes X parameter. Since the Stokes X 
parameter corresponds to intensity, all physically realizable sources will contribute to the 
output through this term. The map of 1Z\ reflects the general imaging behaviour of the bo- 
lometer. This basic output is modified depending on the state of polarisation of the source 
and the values of lZ 2 -4 at the source’s location. 1Z 2 shows how the response is modified if the 
source is a horizontally polarised, IZ 4 if it is diagonally polarised, and IZ 4 if it is circularly 
polarised. 

When p/X » 1, 1Z\ indicates that the device is largely only responsive to the intensity of 
the source over the region bounded by the outline of the re-imaged absorber. lZ\s show that 
the response is slightly polarisation sensitive, reaching 2% of the intensity response at some 
points. However, the polarisation-sensitive regions appear only near the edges of the image. 
Once p/X < 1 the optical behaviour changes significantly. The bolometer becomes responsive 
to the intensity of the field over a circular region that is much larger then the geometrical 
image of the absorber. There is also significant polarisation dependence, which in the case of 
p/X — 0.25 is 10% of the peak intensity response. Notice that the bolometer is insensitive to 
the V component for all values of p/X, which follows from symmetry considerations. Overall, 
these plots indicate that the behaviour changes from being that of a classical optical detector 
to that of a few-mode radio antenna as the size of the absorbing region is reduced below a 
wavelength. 

The increase in polarisation response with decrease in p/X suggests that the state of coher- 
ence of the reception pattern increases, and the number of modes decreases. Figures 4(a)- (c) 
show how the mode coefficients change as p/X is varied. In (a) and (b), the coefficients have 
been scaled by l/(ap 2 ) to remove the dependence on absolute size. In (c) they have been 
scaled by 1/ In to emphasise the relative contribution of each mode to the output. It can 
be seen in (a) and (b) that the number of optical modes to which the bolometer is sensitive 
decreases with p/X. In the limit p « A, the bolometer is sensitive, effectively, to only three 
modes: two of which are degenerate (Figure (c)). In (a) and (b) it can also be seen that for 
p/X < 1 the size of the largest coefficients in each spectrum decreases as p/X decreases. In 
general it was found that, for an ordering based on the magnitude of the coefficients, the 
modes of a given order tend to have same spatial form for all values of p/X. For the modes 
examined, the non-zero part of the field is always confined to a A x A region of the source 
plane centered on the re-imaged pixel. The size of this area relates to the coherence length 
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enforced on the response function on the sky by propagation. 

Figure 5 shows the Stokes response, when p/X — 0.25, broken down into the contributions 
from the three most significant modes. Very similar behaviour is seen for p/X — 0.5. It is 
clear that the polarisation response is dominated by the first two modes. Combining the 
first- and second-mode contributions gives the petal structure seen in Q and U, which is 
enhanced by the addition of the third mode. These two modes are degenerate and therefore 
represent orthogonal states of polarisation of the same spatial mode. Figure 4(c) suggests 
that the behaviour of any electrically small bolometer, p/X < 1, can be modelled by a few 
modes weighted by the appropriate coefficients, and so it is intriguing to consider the forms 
of these modes. 

7. Far-field illumination 

Consider the situation where a bolometer is illuminated directly in the far field. This ar- 
rangement will allow us to explore the state of coherence of the angular reception pattern. 
The field over the surface of the detector can be written in terms of a superposition of plane 
waves: 

M r t,w)= J d% W (k t , u) exp(*k* • r t ), (32) 

where internal scattering has been ignored, and it has been assumed that the reference 
surface is at z = 0. On substituting (32) into (1) and (2), the power absorbed by the 
bolometer becomes 

p(u) = J d 2 k tl J d? k t2 & (k tl ,k t2 ,uj) ■ ■ W(k n , k t2 , u ) , (33) 

where G(k t i, k* 2 , u) is as defined in (10), and VV(k*i, k t2 , u) is as defined by 

W(k tl ,k t2 ,u;)= (W(k u ,cu)W*(k t2 ,o;)). (34) 

The integral is taken over all possible directions of incidence with respect to the front face of 
the reference plane. It may also be true, however, that the bolometer is sensitive to radiation 
incident on its reverse face, and the detector response function can be extended easily to 
include this possibility. We know that the vector k t is related to a wave’s direction of origin 
Cl through 

kt — —k (xx + yy) • Cl — -kl t ■ Cl, (35) 

in addition, 

d 2 kt - -k 2 ( z • Cl)dfl. (36) 
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(33) can then be written 


p(uj) = / dh\ J dfl 2 ^4g(Oi, 5^2 5 w) ■ ■ B(Oi 


^2, w), 


(37) 


where 


and 


= - k A (z • ili) (z • f2o) = . = 

- ofo [L W{-klt • Oi, ~kl t • O^), 

2(27t) 2 Z 0 

^I e (f2i,f2 2 ,^) = 2(2'K) 2 Z 0 Q(-kl t ■ f2i, -kit ■ 


(38) 

(39) 


The factor 2(27 r) 2 Z 0 has been introduced because then fi 2 , a;) takes on the form of 

a brightness tensor, and *4 e (f2i, f2 2 ,cj) the form of an effective-area tensor. This factor can 
be justified more rigorously by using Poynting’s theorem. When the sky is incoherent and 
unpolarised, it may be shown, as expected, that 


and then 


or 


Z3(f2i, CI 2 , (jS) — — {l — B(Ct \ , o;)5(Qi — f2 2 ) 


= i J*t [ HA AA) ] «(*,„), 

p(uS) — J dfl A e (Cl,u ) B(Cl,u>), 


(40) 

(41) 

(42) 


which is a well-known result in radiometry [14] . 

Crucially, it is now possible to find the natural optical modes of the detector in the far- field 
by solving 

J dfi 2 i e (n i; n 2 , u) •uW(fi 2 , w )- 7 <">u< n >(n 1 ,w). (43) 

Expressing the integral in terms of spherical polar coordinates this expression becomes 



sin 62 A e (0i, ^-, 62 ,^ 2 ) 


•Uh9(0 2 ,0 2 ) = 


7 WU^(0 1 ,0 1 ), 


(44) 


where the a;- dependence has been suppressed for brevity. The integrals are taken over the 
full sky, rather then just the forward hemisphere. In its present form, the equation is not 
amenable to solution by the numerical, modal method discussed previously. The factor sin 0 2 
means the that the kernel is non-symmetric, unlike the underlying function. Instead, consider 
the change of variable u — — cos 9. With this substitution, (43) becomes 



d<t >2 MOM, 0 1 ; 6 (u 2 ), (j> 2 ) • U<">(0(ti2), <h) 


AWuW(0( Ul ),^). 


(45) 
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In this form, the modes, U ^(9(u 2 ), ^2), and mode-coefficients, 7^, can be found by sampling 
at regularly spaced intervals in the u- and ^-ordinate, and finding the eigenvalues and ei- 
genvectors of the matrix AuA(j}A e {0{ui), </>i; 9(u 2 ), fa)- Some care must be taken with the 
values of A (f> and A u. In the case of the ^ordinate, the integral is taken around a closed 
loop and as a result the A is constant. 


8. Simulations II 

The above method was used to calculate the natural modes of the angular reception pattern 
of the bolometer described in Section 4. Substituting k t2 , to) for the model pixel into 

(39), we find the angular response of the model pixel is given by 


A e (&i, & 2 ;to) = 2ap 2 Z 0 l sine ^(^1 — D 2 ) • x 

-A 


p 

x sine L (n, - ri 2 ) • yj , 


(46) 


where we assume the response in the reverse direction is the same as in the forward direction. 
Figure 6 is the equivalent of 4(c) for the far- field mode coefficients. As in the case of re-imaged 
pixels, the number of significant modes decreases rapidly as p/X falls below unity. The number 
of modes in this limit tends to three, which are degenerate in this case. We also calculated, 
although they are not shown here, the angular forms of the first four modes when p/X — 0.25. 
The modes closely resemble electric and magnetic dipoles having different orientations. To 
establish equivalence, we attempted to fit electric- and magnetic-dipole fields to each of the 
first four modes. The first three modes showed excellent agreement with fields for x-, y- and 
z-orientated electric dipoles, and the fourth mode to a z-directed magnetic dipole. 


9. Dipole representations 

At the end of the previous section we observed that when p/X — 0.25 the far field behaviour 
is dominated by a set of dipole- like modes. This suggests that for small p/X, (46) should be 
well described by 


JUft„si 2 ,u,) ~ ]T £ D„,(n 1 )D;,(0 2 ) 


(47) 


where (D n (D)} is the set of spherical vector harmonics corresponding to electric- and 
magnetic-dipoles orientated in the x, y and z directions. The matrix elements are given 

by 

AnnM = J dtii f dtl 2 D^(Di) • n 2 ,u) • D„(fi 2 ). (48) 

If A mn is diagonal, it does not automatically follow that the optical modes correspond 
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to dipole fields. A mode could, for example, consist of a dipole field in combination with a 
higher order spherical harmonic. A measure of how good the approximation is at a particular 
wavelength can be obtained by calculating the ratio A nn / where we have, using 

(46). 

^^ 7 m — [ d£l X • • A e (f2, Cl,u) — 16nap 2 Z 0 . (49) 

m 

If the dipoles correspond to the full set of optical modes, then ^2 n A nn — Another 

useful quantity is A nn j s fa Jm 7 m , which gives the relative responsivity of a bolometer in a 
particular dipole mode. 

The challenge is to find the matrix elements A mn (uj). The reception patterns associated 
with differently orientated magnetic dipoles are given by 


M a .(0, 0) — N ^sin 00 + cos 9 cos 0 (frj 
M 2/ (0, 0) — N (— cos 0 0 + cos 9 sin 0 <0 j 
M 2 ,(0, 0) — N sin 9 0, 


(50) 


where 0 and 0 are unit vectors, and N is the normalisation factor i/3/ (8it). Using the same 
notation, the fields associated with differently orientated electric dipoles are 

E x (9, 0) — N (cos 9 cos 0 9 — sin 0 

E y (0, 0) = AT ^cos0sin00 + cos00^ (51) 

E z (9, 0) — N s\u 9 9. 

When using (50) and (51) in (48) it is beneficial to introduce the operator 

: / A A A ^ \ / A A A^\ 

T — (0101 + 0101 J • [0202 + 0202 J 

= [cos 9\ cos 92 cos(0i — fa) + sin 9\ sin 0 2 ] 0102 + cos 9\ sin(0i — fa) 9\fa (^2) 
+ COS 9 2 sin (02 - fa) 01^2 + COS (0i - fa) fa fa- 


Strictly speaking, each of the brackets in the first line of (52) also contains a term of the 
form rr ; however, since neither the magnetic dipole fields or the electric dipole fields have 
a component in the radial direction this term can be dropped without effect. Replacing the 
idem factor in (46) with T is beneficial because the dipole fields can be represented as two- 

_ ^ /s 

element vectors. The conversion between [0i,0 x ] and [02, fa] is then handled in a by T, 
which can be represented by a 2 x 2 matrix. 
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In spherical coordinates, (48) becomes 

/* 7T 7*77 /» 27T p2-k _ 

A mn(u) = dO 1 do 2 /#l /# 2 Sm 01 sin 02 D^(6>i , 0i) • AOA 0i; # 2 , A • D n (0 2 , A)- 
Jo Jo Jo Jo 

(53) 

This expression is difficult to solve analytically, and therefore the integrals were evaluated 
numerically. To check the software, the limiting case p — > 0 was considered. As p — > 0, from 
(48) and (46) it can be shown that 

A (lu) 1 r* f 2n f 2 * 

16? 7ap 2 Z 0 ~ ' 8 tt J Mi J^ 2 * J d( ^ 2 sin6> i sin0 2 D^(6>i,0i) • D n (6> 2 ,0 2 )- (54) 

The integrals that arise in this case can be evaluated analytically. It is found that only the 
diagonal elements associated with the electric dipole modes are non-zero, each being equal 
to 1/3. These analytical limits were used to check the software. 

10. Simulations III 

Figure 7 shows how the diagonal elements, A nn , vary with p/X. Instead of plotting the A nn 
directly against p/X, the scaled values A nn /l67iap 2 Z 0 are shown. There are two advantages 
with this normalisation: The first is that the curves become independent of the bolometer’s 
sensitivity a and its absolute size p; the coefficients of a particular device can be found 
by reading off the graph and rescaling. The second advantage is that the curves provide a 
measure of the relative responsivity to each dipole mode. The numerical error in the values is 
of the order of ±1 x 10~ 4 . It can be seen from Figure 7 that the numerical values of A nn in the 
limit p/X —> 0 are in agreement with the analytic values discussed at the end of the previous 
section. In addition, the values derived through this method are in agreement with those 
that were found when carrying out the full simulations in Section 8. For non-zero p/X, it was 
found that the off-diagonal elements, A mn , are zero (within numerical error) for all values 
of p/X < 1.5. This observation indicates that the reception pattern corresponds, to a good 
order, to an incoherent superposition of electric and magnetic dipoles. The values shown in 
Figure 7 are proportional to the responsivity of the detector in each of the dipole fields. It 
can be seen that the bolometer is equally responsive to power in the x- and y-orientated 
electric dipole fields, and similarly in the x- and y-orientated magnetic dipole fields. This 
degeneracy is expected from the symmetry of the absorber. As p/X falls, the responsivity 
to the electric dipole fields increases, and tends to a limit at p/X — 0. The responsivity to 
the magnetic dipole fields increases, peaks atp/A«0.75, and then decreases to zero, as the 
area tends to zero. The ratio Y/ n Am/ Y2 m 7™ a particular value of p/X may be calculated 
from Figure 7 by summing the contributions from each line. Over the range 0 < p/X < 0.5, 
Aa -Am/ A m 7 m is close to unity. In this region, the dipole representation therefore provides 
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a reasonably complete description of the bolometer’s behaviour. 


11. Thin resistive films 


To this point we have considered volumetric absorbers, but the methods described can be 
extended to thin resistive films. ‘Thin’ is used here to mean the thickness is very much 
smaller then the skin depth, so currents are confined to flow in the plane of the film. In 
operation, an incident field induces a surface current density J s on the film. Ohmic heating 
by these currents then results in the film extracting power from the incident field. If the film 
has surface resistance R s , and is square with side length p, the absorbed power is 


p(u) 


1 

2 



P 

2 


dy R s \J a (x,y)\ 


P 

2 


2 


(55) 


The surface current density results from the total field, and therefore contains contributions 
from both the incident field and the field scattered by the currents. Strictly, the scattered 
field and surface current density are related via the Green’s dyadic for the space, so a 
rigorous calculation of the induced current requires solving an integral equation. However, 
by assuming the induced surface current is a uniform, infinite, sheet an approximate linear 
relationship between the two quantities may be found. This approach ignores the effects of 
the edges of the film and the direction of incidence of the radiation, but is sufficient for a 
first-order model. It can be shown [15] that the field generated by an infinite sheet current 
lying in the plane z — 0 is 

E s(z) = -yJ*exp(ifc|*|). (56) 

Taking (56) as the scattered field and using the constitutive equation Z a J s — E (z) + E s (z), 
we find that the induced current can be written in terms of the incident field E(z) as 


J s(x,y) 


1 

(■ Rs + Zq/2) 


(: I - zz) 


•E (x,y). 


(57) 


The dyadic operator picks out the component of the incident field tangential to the surface, 
which is necessary since currents can only flow in the plane of the film. (57) shows that the 
physical effect of the scattering is to increase the apparent surface impedance of the film. 
Substituting (57) into (55), we obtain 


1 


pH =o 


R s 


2 (R s + Zq/2)- 


P P P P 

/ dx i / dy i / dx 2 / dy 2 
J-2 J-2 J-2 J-2 


(l - zz) - (l - zz) 5(xx - x 2 )5{yi - y 2 ) f • • F,(x 1 ,y 1 )E*(x 2 ,y 2 ), 


(58) 
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which, when compared with (2), suggests a planar response function of the form 


R s 


V{v n -r n-u) = { HR s +z 0/ 2 ? 


(l - Z1Z1 ) • (x - z 2 z 2 ) 6 (r u - r t2 ) r a and r t2 on film 
0 otherwise 


(59) 

With the exception of the idem factor, (59) is identical to (20) with a — R S /2(R S + Zq/2 ) 2 , 
which is the power mismatch factor. This mismatch factor could also have been calculated 
on the basis of transmission line theory. The dyadic operator in (59) is position independent, 
and so can be taken out of the integral in (10). The angular response function of the resistive 
film is therefore simply (46) with the X replaced by {x — z x z ^ • (x — z 2 z 2 j . In a spherical 
coordinate system, this dyadic constant becomes 



ZiZi 


(x — z 2 z 2 ^ - cos 6\ cos 9 2 cos(0i — 0 2 ) 0\0 2 + cos 9\ sin(0i — 0 2 ) 0i4> 2 
+ cos 9 2 sin (0 2 - 0i) <^j0 2 + cos (0i - 0 2 ) ^ x 4> 2 , 


(60) 


where, as in (52), we have suppressed the parts involving ?, as dipole fields have no component 
in the radial direction. As <0 has no component in the z direction, (60) differs from (52) only 
in the term in 0\0 2 . Consequently, a simple modification is all that is required in order to 
model thin-film devices. 


12. Simulations IV 

Figure 8 shows how the diagonal elements of A mn / Y! n ln of the thin-film bolometer vary 
with p/X. For the resistive film bolometer J2nXn — (I6np 2 Z 0 R S ) /[3(R S + Z 0 /2) 2 ]. The off- 
diagonal elements have not been shown, as they are all zero within numerical error. Unlike 
the volumetric bolometer, the thin-film bolometer is unresponsive to power in the x- and y- 
orientated magnetic dipole fields and the ^-orientated electric dipole field for all values of p/X. 
This behaviour can be explained by considering the current distributions. The ^-orientated 
electric dipole field couples to a sheet of current in the z-direction, while the x- and y- 
orientated magnetic dipole fields couple to loops of current in the (x, z ) and (y, z ) planes. All 
of these distributions require a current normal to the surface of the film, and therefore cannot 
be supported. The modes that are supported — the x- and y-orientated electric dipoles and 
the z-oriented magnetic dipole — couple to currents in the plane of the film. For the modes 
that are supported, the behaviour with changing p/X is the same as for the volumetric case. 
Similarly, Am/ obtained by summing the lines at a particular value of p/X, is 

greater than 0.8 for p/X < 0.5. Again this indicates the dipole representation provides a good 
description of the bolometer’s far field optical behaviour for small p/X. 
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13. Conclusion 


We have considered Withington’s [6] and Saklatvala’s [7] theoretical model of the optical be- 
haviour of free-space power detectors. In the first part of the paper, we presented an efficient 
numerical technique for simulating the behaviour of ideal bolometric imaging systems, and 
demonstrated how the intensity and polarisation response change as the size of the absorb- 
ing region is reduced below a wavelength: the reception pattern changes, effectively, from 
being that of a highly multimode classical optical detector to that of a few-mode radio an- 
tenna. The changing form results largely from the limitations imposed by the free-space path 
between the source and detector, and therefore the observed behaviour should be character- 
istic of all volumetric absorbers. In the second part of the paper, we extended Withington 
and Saklatvala’s theory, and considered the angular response of naked pixels. It was found, 
for volumetric absorbers having a side length of less than a wavelength, six modes dominate 
the intensity and polarisation response: three orthogonal electrical dipoles and three ortho- 
gonal magnetic dipoles. As the size of the absorber changes, the relative contributions of the 
modes change. The scheme can be modified easily to cover thin-film absorbers. In this case, 
current can only be induced in the plane of the film, and the response of small bolometers 
is dominated by two electric and one magnetic dipoles. By using a striped, rather than con- 
tinuous, resistive film, one would expect to be able to reduce the response to a single dipole 
mode. Not only does this work demonstrate that the optical behaviour of any bolometer 
can be modelled straightforwardly, it also suggests a particularly simple way of modelling 
the behaviour of electrically small bolometers. The ability to simulate easily the imaging, 
polarisation, and straylight characteristics of small pixels will be of considerable value when 
designing imaging arrays. 
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f = focal length 



F s 

Fig. 1. An ideal imaging system with a free-space bolometer. 
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Fig. 2. Figure (a) shows the sampling points in the ^-domain, while figure 
(b) shows the corresponding sampling points in the spatial domain. In each 
case the total number of points is the same. The grey area in (a) indicates the 
region of the fc-plane where the modes are non-zero. 
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Fig. 3. The Stokes parameters of the optical system shown in Figure 1 for 
values of p/A in the range 0.25-4.0. In each plot, the positive x-direction is 
horizontal from left to right, and the positive y-direction vertical from bottom 
to top. The small black square is the re- imaged pixel. For each value of p/A, all 
four plots are normalised to the peak flux in the corresponding X plot. Each 
set of U , Q and V plots use the same colour scale, but the scale’s range changes 
with p/A. 
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Fig. 4. Spectra of mode coefficients for different values of p/X. In (a) and (b) 
the coefficients, 7 have been scaled by l/(aip 2 ); in (c) the coefficients have 
been scaled by 1 / Yhln to indicate relative contribution. 
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Fig. 5. The contributions from the three coherent-modes with the largest 
eigenvalues to the Stokes parameters of the reception pattern when p/ A=0.25. 
The small black square is the re-imaged pixel. All plots are normalised to the 
peak flux in the Stokes X of the full field. The Stokes parameters are plotted 
against a single colour scale to emphasise their relative contributions. 



Fig. 6. The normalised coefficients, 7n/ An, of the far-field optical modes. 
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Fig. 7. Diagonal elements of A mn as a function of p/X for a volumetric ab- 
sorber. The ratio Yl n Ann/ 7m can ^e calculated for a particular value of 
p/X by summing the contributions from each line. This ratio provides a meas- 
ure of how well the dipole model describes the behaviour of the absorber, and 
is equal to unity when the description is complete. 



Fig. 8. Diagonal elements of A mn as a function of p/X for the resistive film 
absorber 
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